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I.  INTRODUCTION 


References  1  through  6  have  demonstrated  that  the  one-dimensional  (1-D)  plane  impact 
shock  sensitivity  of  several  important  energetic  polymeric  chemical  CHNO  compounds  can  be 
correlated  to  how  much  heat  energy,  A  (v.e.)iR,  they  can  absorb  before  some  type  of  reaction 
(melt,  phase  transition,  deflagration,  or  detonation)  occurs.  The  amount  of  internal  energy  that 
can  be  soaked  up  is  directly  related  to  the  material’s  specific  heat  (Cp)  magnitude  and  variation 
with  temperature  (T).  See  Section  II. A.  for  more  details  about  the  A  (v.e.)xR  concept  that  relates 
Cp  heat  absorption  to  impact  shock  sensitivity. 

The  Cp  (T)  of  CHNO  energetic  polymers  is  a  very  important  property  from  both  heat  and 
impact  shock  sensitivity  view  points.  As  such,  a  simple  Cp  (T)  estimation/prediction  procedure 
is  also  desirable  and  very  important  when  experimental  data  may  not  be  available. 

One  attempt  at  a  Cp  (T)  estimation/prediction  scheme  is  documented  in  Reference  7  and 
briefly  described  in  Section  II. B.  of  this  report.  This  is  the  so  called  Nominal/Generic  (N/G)  Cp 
per  average  atom  concept. 

An  N/G  Cp  was  proposed  based  on  the  observation  that  the  CP’s  per  average  atom  for 
CHNO  explosives  at  a  given  temperature  did  not  deviate  very  much  from  each  other.  So  in  lieu 
of  experimental  Cp  data  for  a  CHNO  explosive  compound,  the  N/G  Cp  could  be  utilized  to  make 
a  Cp  prediction. 

A  successful  application  of  this  Cp  estimation  procedure  is  documented  in  Section  III  of  this 
report  for  the  FOX-7  CHNO  explosive  compound.  This  application  is  called  “successful” 
because  the  predicted  Cp  agreed  remarkably  well  with  Cp  results  computed  by  a  much  more 
complex  and  sophisticated  analysis. 
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II.  BACKGROUND  INFORMATION 


A.  The  A  (v.e.)xR  Concept 

Essentially,  the  area  A  (v.e.)iR  under  the  Cp  versus  Temperature  (T)  plots  between 
temperature  limits  (Experimental  Test  Temperature  (Texp)  and  Reactive  Temperature  (Tr))  is  a 
measure  of  how  much  atomic  vibratory  energy  explosives  can  absorb  before  a  reaction  occurs. 
The  reaction  may  be  melting,  phase  change,  decomposition,  burning,  or  even  detonation.  Thus, 
to  a  good  approximation,  it  could  be  expected  that  if  A(v.e)iR  amount  of  energy  is  suddenly 
added  via  impact  shock  loading,  then  a  reaction  may  occur. 

This  A  (v.e.)xR  concept,  that  impact  shock  sensitivity  or  shock  induced  reactivity  of 
energetic  materials  could  be  related  to  their  specific  heat  (Cp)  variation  with  temperature,  was 
demonstrated  in  References  1  and  2  for  RDX,  TETRYL,  PETN,  TNT,  and  TATB,  which  are 
basic  secondary  reactive  compounds. 

In  References  3  and  4,  the  A  (v.e.)xR  ideas  were  demonstrated  for  HMX  and  HNS 
which  are  also  important  basic  secondary  explosive  compounds.  The  impact  shock  response  of 
these  seven  compounds  ranges  from  very  insensitive  to  highly  sensitive.  Most  of  these  seven 
basic  energetic  compounds  have  been  the  main  ingredient  of  useful  explosive  mixtures. 

One  such  mixture  is  the  plastic  bonded  explosive  designated  as  PBX-9502  that  is  95 
percent  TATB  and  5  percent  KEL-F80.  PBX-9502  has  been  rather  extensively  tested  via  1-D 
shock  loading  at  various  temperatures,  and  its  thermal  characteristics  have  also  been 
experimentally  explored.  Consequently,  with  this  much  information  available,  A  (v.e.)xR 
concept  computations  were  made  for  PBX-9502.  The  exploratory  comparative  results  for  this 
important  energetic  material  were  highly  affirmative  and  are  documented  in  References  5  and  6. 

General  details  of  the  exploratory  computation  and  experimental  data  comparisons 
involved  in  a  general  A  (v.e.)xR  assessment  are  contained  in  the  following  paragraphs. 

The  thermal  atomic  vibratory  energy  increment,  A  (v.e.)xR,  is  related  to  macroscopic 
critical  particle  (or  mass)  velocities  (  Upcr  )  and  impact  shock  pressures  via  Equations  (1) 
through  (6). 

For  certain  explosives,  a  good  estimate  of  the  critical  particle  velocity,  Upcr  ,  where  a 
reaction  (or  detonation)  occurs  is: 


T  T  _  A(v.e.)iR 

U  PCRl  —  ^  - 


mAv 


In  some  circumstances,  a  better  estimate  of  the  critical  particle  velocity  is: 


Upc 


2A(v.e.)r 


mAv 


V2  Upc: 


(1) 


(2) 
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where: 


A(v.e.)TR 

CP 

ITLav 

Texp 

Tr 

UPCR2 

Upcri 

et 

e; 


-  C  c>dT  <3> 

=  Thermal  vibratory  energy  per  atom  between  TExp  and  Tr, 

Gram  (Cm/Sec)2. 

=  Specific  heat  per  atom  as  a  function  of  temperature. 

=  Average  mass  of  an  atom  in  the  material,  Grams  (Appendix  A). 

=  Temperature  at  which  experimental  impact  shock  tests  are 
conducted.  This  is  normally  room  temperature  (RT~  300  °K) 
but  can  (and  should)  be  done  at  higher  and  lower  temperatures. 

=  Temperature  at  which  some  thermally  induced  reaction  occurs 
(decomposition,  melting,  phase  change,  detonation,  etc.). 

=  Particle  velocity,  Up,  such  that  the  shock  induced  internal 
energy  (e,)  is  equal  to  A(v.e,)iR. 

=  Particle  velocity,  UP,  such  that  the  total  shock  induced  energy 
(et)  (kinetic  plus  internal)  is  equal  to  A(v.e.)iR. 

=  mAvU  p  =  total  shock  energy  per  average  atom. 

=  ^U'  =  ek  =  internal  or  kinetic  energy  of  the  shocked 

material  per  average  atom. 


When  good  comparisons  can  be  made,  it  has  been  generally  found  that  experimental 
detonation  particle  (mass)  velocity  threshold  values  (  Up™  )  will  be  somewhere  between  Upcri 

and  Up  CR2  . 


Certain  explosives,  when  heated  to  higher  and  higher  temperatures,  melt  before  they 
explode  (RDX  and  TNT,  for  example).  This  melting  will  require  that  the  heat  of  fusion  (AHp)  be 
absorbed  by  the  material  at  T=TMelt  conditions  before  the  temperature  will  increase  [9,  10]. 
Consequently,  if  TExp  is  less  than  Tmelt,  then  the  total  heat  absorbed  from  T  =  Texp  to  T  =  Texpl 
=  Tr  is: 


c  Tmelt  /•Texpl 

A(v.e.)TR  =  CpdT  +  A  HF  +  [  CPdT.  (4) 

Therefore,  for  solid  energetic  materials  which  melt  prior  to  explosion,  then  A  (v.c.)ir, 
as  defined  by  Equation  (4)  is  employed  in  Equations  (1)  and  (2)  to  compute  Upcri  and  UPCR2, 
respectively.  Note  that  melting  is  just  one  example  of  a  phase  transfonnation  which  may  require 
an  enthalpy  increment  (AHj)  to  be  activated.  For  example,  HMX  can  exist  in  different  solid 
polymorphic  forms.  At  a  certain  temperature,  Tt,  one  form  may  change  to  another  fonn  if  the 
heat  energy  of  transformation  (A  HT)  is  supplied.  Therefore,  AHT  should  be  added  to  Equations 
(3)  and  (4)  if  Tr  is  greater  than  T-p. 
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Note  that  A  (v.e.)iR  as  defined  by  Equations  (3)  and  (4)  is  actually  an  enthalpy 
increment  (AH).  However,  it  was  shown  via  numerical  examples  in  Appendix  B  of  Reference  1 
that,  under  the  experimental  Cp  acquisition  conditions,  the  pressure  times  volume  tenns  were 

pTR  pTR 

minute  compared  to  the  Cp  integral,  CpdT.  Thus,  CpdT  is  essentially  all  of  the 

JTexp  JTexp 

internal  energy  difference  caused  by  thermal  stimulation  during  standard  tests  at  atmospheric 
pressure  to  determine  the  specific  heat  characteristics. 

Once  A  (v.e.)iR  and  Upcr  values  are  computed,  the  corresponding  shock  velocity 
(Uscr)  is  ascertained  from  experimental  data  for  Us  as  a  function  of  the  particle  velocity,  Up. 
The  experimental  relationship  is  usually  linear  and  written  empirically  as: 

Us  =  Co  +  SUP  (5) 

When  Up  =  Upcr  and  Us  =  Uscr  are  determined,  the  shock  pressure  is  computed  from 
the  following  well  known  relation: 

Ps  =  fx  Us  Up,  (6) 

Where  p„  =  Material  density  (grams/cm  ). 

Then  Upcr  ,  Uscr  ,  and  Pscr  may  be  compared  to  experimental  shock-induced  reaction 
threshold  information  to  check  the  validity  of  the  above  A  (v.e.)iR  theory  to  denote  reactive 
conditions  under  impact  shock  stimuli.  The  numerical  computations  involved  in  a  A  (v.e.)xR 
assessment  are  straightforward  and  may  be  perfonned  with  a  hand-held  calculator. 

It  must  be  emphasized  that  any  possible  effect  of  pressure  on  Cp  is  not  taken  into 
account  in  the  present  analysis.  The  basic  idea  is  that  if  a  quantity  of  thennal  vibration  energy, 
A  (v.e.)iR,  under  quiescent  conditions  is  able  to  create  a  reaction,  then  the  same  amount  of 
energy  added  by  an  impact  shock  (e;  or  et)  should  also  cause  some  type  of  reaction.  The  shock- 
induced  reaction  may  not  be  the  same  type  as  the  temperature  induced  reaction,  but  will 
nevertheless,  be  a  reaction  of  some  kind.  It  may  be  less  or  more  severe  than  the  thermally 
induced  reaction. 

The  Cp  unit  of  calorie/(atom  °K)  was  employed  in  plots  of  Cp  versus  T  information, 

23 

which  are  shown  in  this  report.  This  is  because  the  Boltzman  constant,  kB=0.33  x  10'“ 
calories/(atom  °K)  and  the  maximum  Cp  at  high  temperatures  for  many  materials  is  3  1<b  ~  1.0  x 
10'23  calories/(atom  °K).  This  is  a  good  mnemonic  reference  level  for  comparison  purposes.  It 
was  noted  in  Reference  1  that  the  average  Cp  per  atom  for  most  polymers  never  reaches  the  3  ks 
level  before  a  reaction  (phase  change,  melting,  glass-to-rubber  transition,  or  even  detonation) 
occurs. 


Actually,  Cp  for  some  atoms,  or  combinations  of  atoms,  probably  reaches  the  3  ke 
level  and  causes  a  reaction  at  some  Tr.  But  Cp  for  a  large  number  of  atoms  remains  much  less 
than  2  ke.  Thus,  a  large  amount  of  the  possible  thermal  vibratory  energy  is  never  activated  and 
the  average  Cp  per  atom  remains  relatively  low  [14].  In  many  cases,  important  temperature 
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induced  reactions  occur  near  the  average  Cp  *  2  1<b  level  at  moderate  temperatures  (400  to 
600  °K). 

B.  The  N/G  Cp  Concept 

It  was  first  documented  in  Reference  1  that  the  Cp  [Cal/(atom  °K)]  magnitudes  (at  a 
given  T)  for  five  solid  explosive  compounds  (RDX,  TETRYL,  PETN,  TNT,  and  TATB)  did  not 
differ  very  much  from  each  other.  In  Reference  3,  it  was  demonstrated  that  the  Cp  for  HNS  was 
very  close  to  that  for  TATB,  and  the  Cp  for  HMX  was  somewhat  less  than  the  TATB  Cp  at  the 
higher  temperatures.  Liquid  TNT  has  a  larger  CP  (at  a  given  T)  than  these  solid  energetic 
compounds.  These  statements  are  corroborated  by  the  experimental  Cp  information  exhibited  in 
Figures  1  through  4  of  Reference  7. 

So  with  two  exceptions  (melted  TNT  and  5  -  HMX),  the  Cp  per  atom  of  five 
important  secondary  energetic  compounds  all  had  very  similar  magnitudes  near  that  for  TATB. 
HNS  also  has  the  same  magnitude  and  trend  (non-linear  variation)  with  temperatures  as  TATB. 
The  other  explosives  (RDX,  TNT,  TETRYL  and  PETN)  Cp  have  a  linear  variation  (a  +  b  T)  over 
most  of  their  temperature  range. 

The  amount  of  heat  energy  these  compounds  and  mixtures  can  absorb  varies 
considerately.  TATB  and  PBX-9502  can  soak  up  more  heat  energy,  A  (v.e.)iR,  than  the  others 
by  a  considerable  margin.  That  is,  their  reactive  temperatures,  TR,  were  much  greater  than  the 
other  compounds  (more  than  100  °K  for  TNT,  HNS  and  HMX  and  over  200  °K  for  TETRYL, 
PETN,  and  RDX). 

So  based  on  the  above  remarks  and  similar  remarks  in  References  1,3,  and  5,  a 
nominal  Cpper  average  atom  was  proposed  and  defined  as  shown  in  Figure  2  and  listed  in 
Table  1.  These  nominal  Cp  magnitudes  at  Room  Temperature  (RT)  and  above  are  very  similar 
to  those  of  TATB,  particularly  at  the  high  temperatures.  Below  RT,  the  proposed  N/G  Cp 
magnitudes  are  very  close  to  (or  equal  to)  the  RDX  Cp  values  and  TNT  Cp  values  near  absolute 
zero. 


In  Reference  7,  two  examples  (TNT  and  HMX)  were  selected  for  comparative 
computations  where  Upcr  andPscR  were  calculated  via  the  experimental  Cpand  the  N/G  Cp.  TNT 
has  a  high  Cp  relative  to  the  N/G  Cp  and  HMX  has  a  low  Cp  relative  to  the  N/G  Cp.  Two  cases 
were  considered  for  TNT  where  the  large  liquid  CP  was  included  in  Case  1 . 

Consequently,  excluding  Case  1  for  TNT  where  the  large  liquid  Cp  had  to  be  included 
in  the  Upcr  (EXP)  and  Pscr  (EXP)  computations,  then  for  TNT  (Case  2),  PBX-9502,  and  HMX, 
the  maximum  percentage  differences  between  the  (EXP)  CP  and  (N/G)  CP  results  were  bounded 
by: 


AUpcr 

=  UMN/G)-  UMEXP)  x  100.0  <  I  2.5%  1 

(9) 

Upcr(EXP) 

Upcr(EXP) 

APscr 

=  P-(N/G)-PMEXP)X1000 

(10) 

Pscr(EXP) 

Pscr(EXP) 
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These  small  differences  provide  considerable  credibility  for  the  proposed  N/G  Cpper 
average  atom  concept  for  most  CHNO  energetic  materials. 

III.  EXAMPLE  OF  CP  PREDICTION  FOR  FOX-7  VIA  THE  N/G  CONCEPT 

The  relatively  new  CHNO  explosive  compound,  FOX-7  (1,  1 -Diamino-2,  2- 
Dinitroethylene),  has  generated  a  considerable  amount  of  interest  among  explosive  investigators 
within  the  last  decade.  References  8  through  1 1  are  a  good  sampling  of  the  experimental  and 
theoretical  efforts  devoted  to  this  energetic  material. 

However,  there  has  been  no  experimental  detennination  of  FOX-7’ s  specific  heat,  but  there 
is  a  rather  sophisticated  and  complex  theoretical  prediction  of  Cp  as  a  function  of  T  given  in 
Reference  8.  As  such,  this  fits  our  search  criteria  for  a  good  CHNO  energetic  compound  to 
apply  the  N/G  CP  concept  for  predicting  the  CP  magnitudes. 

Thus,  FOX-7  is  a  good  choice  for  a  Cp  prediction  example  Via  the  N/G  Cp  concept  because: 

1.  No  experimental  Cp  (T)  data  for  FOX  -7  were  available  and,  consequently,  could  not 
have  influenced  the  proposed  N/G  Cp  magnitudes. 

2.  There  was,  however,  a  theoretical  computation  of  Cp  (T)  for  FOX-7  that  was 
available  for  comparison  with  Cp  (T)  from  any  other  predictive  scheme  or  with 
experimental  data. 

In  Appendix  A,  the  Molecular  Weight  (MW)  and  mAv  (mass  of  an  average  atom)  for  FOX-7 
are  computed.  These  values  are: 

MW  =  148.086  grams/mole 

mAv  =  1.756239  (10  23)  grams/atom 

Also  for  C2  H4  N4  04  =  Q  Hj  Nk  Oi 

q  =  i+j  +  k  +  l=  14.0  atoms/molecule 

From  various  sources,  Avogadro’s  number  is: 

Nav  =  6.02252  (10+23)=  molecules/MW 

=  Number  of  molecules  per  mole  (MW)  in  a  chemical 

compound.  Also  NAV  is  the  number  of  atoms  in  an  Atomic 
Weight  (AW)  of  an  element. 

As  shown  in  Appendix  A,  the  number  of  atoms  per  mole  (MW)  of  FOX-7  can  be  computed 
by  two  different  ways.  These  are: 
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MW 


148.086 


1.  NAPMWi 


2.  NAPMW2 


niAv 


Grams/MW 
Grams/Atom  1.756329(10  23 ) 


84.315638  (10+23)  Atoms/MW 


,T  Atoms  Molecules 

q  NAv  =  -  *  - 

Molecule  MW 

(14.0)  (6.02252  x  10+23) 

84.315280  (10+23) 

MW 


These  results  for  NAPMW  agree  to  three  decimal  places  and  so 


NAPMW  =  84.315  (10+23)  Atoms/MW 


is  used  in  the  following  computations  to  predict  the  FOX-7  Cp  from  the  N/G  Cp.  That  is,  for  a 
given  temperature  (T): 


Also: 


FOX-7  CP 


FOX-7  CP 


(NAPMW)  Cpn/g 


/  Atom  ) 

v  MW  J 


4.184  CP 


Cal.  ) 

(  Cal.  3 

Atom  -  °K, 

[mw-°k  J 

Cal.  "j 

r  Joules  ^ 

MW-°kJ 

Imw-°k  J 

These  Cp  prediction  computations  are  listed  in  Table  2  for  a  wide  range  of  temperatures 
(100  to  650  °K).  The  results,  Cp  [Joules/  (MW-°K)],  are  plotted  versus  T  (°K)  in  Figure  1. 


Table  3  lists  the  FOX-7  theoretical  Cp  [Joules/  (MW-°K)]  results  from  Reference  8,  (Fig.  6). 
This  infonnation  is  also  plotted  in  Figure  1  for  comparative  purposes.  The  Cp  magnitudes  from 
these  two  predictive  schemes  are  in  remarkably  close  agreement. 


Table  3  also  contains  the  conversion  of  the  Reference  8  FOX-7  Cp  in  [Cal/  (Atom-°K)]. 
These  results  are  plotted  in  Figure  2  along  with  the  proposed  N/G  Cp  that  is  in  the  same  units. 

As  expected  from  the  Figure  1  comparison,  the  FOX-7  CP  and  the  N/G  CP,  on  a  per  average  atom 
basis,  agree  exceptionally  well. 
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IV.  DISCUSSION 


As  shown  in  Figure  1,  the  FOX-7  Cp  (T)  magnitudes,  predicted  from  the  proposed  N/G  Cp 
concept  in  the  present  report  and  the  theoretical  Cp  results  in  Reference  8,  compare  exceptionally 
well.  Table  4  lists  the  magnitude  and  the  percentage  differences.  The  percentage  difference  is 
within  6.0  percent  for  100  °K  <T<  400  °K. 

The  largest  negative  difference  (-5.9  percent)  occurs  at  100  °K  and  the  greatest  positive 
difference  (+3.9  percent)  occurs  at  400  °K.  So  the  predictive  trends  (or  slopes)  with  temperature 
are  slightly  different  but  the  Cp  magnitudes  are  quite  comparable.  This  is  remarkable, 
considering  the  differences  in  sophistication  and  complexity  between  the  predictive  methods. 

V.  RECOMMENDATIONS 

As  a  result  of  the  present  analysis  and  comparison  of  predicted  Cp  (T)  results  for  the  FOX-7 
CHNO  explosive  compound,  the  following  recommendations  are  made: 

1 .  Experimental  Cp  data  for  FOX-7  should  be  acquired: 

a.  To  compare  with  both  Cp  predictions. 

b.  To  have  the  authoritative  Cp  for  use  in  the  A  (v.e.)xR  concept  impact  shock 
sensitivity  computations. 

2.  The  temperature  range  should  be  from  100  °K  to  a  practical  upper  limit  from 
experimental  considerations.  See  Reference  1 1  for  documentation  of  early  work 
on  thermal  sensitivity. 

3.  Plane  shock  impact  experiments  should  be  performed  to  ascertain  a  lower 
threshold  detonation  limit  for  UP  or  Ps.  This  should  be  done  for  more  than  one 
temperature  (TExp).  Reference  5  documents  a  A  (v.e.)iR  concept  shock  sensitivity 
analysis  for  PBX-9502  at  four  Texp  temperatures  (218,  300,  348  and  525  °K). 

4.  Predict  the  FOX-7  Cp  (T)  variation  using  Satoh’s  scheme  for  submolar  group 
additive  contributions  to  the  total  Cp  of  a  polymer  compound.  See  References  12 
and  13  for  further  infonnation  on  this  additive  property  as  applied  to  inert 
polymers. 

With  regard  to  recommendations  1,  2  and  3,  this  experimental  data  may  be  available  for 
FOX-7,  but  if  so,  it  is  not  known  to  the  author  of  this  report. 

With  regard  to  recommendation  No.  4,  some  rather  affirmative  preliminary  Cp  predictions 
via  the  Satoh  submolar  group  additive  property  have  been  made  for  TNT  and  RDX.  Similar 
computations  will  be  made  for  FOX-7  and  all  these  Cp  predictive  results  for  energetic  polymer 
compounds  will  be  published  in  a  forthcoming  report. 

The  above  recommendations/applications/comments  concerning  explosives  in  this  report 
should  also  be  valid  for  solid  CHNO  rocket  propellants. 


Table  1 .  The  Generic  Cp  per  Average  Atom  for  CHNO  Explosives 


T 

T 

CP 

Remarks 

°K 

°C 

Cal 

Atom  -  °K 

~ 

0.0 

-273 

0.0000 

t 

7.0 

-266 

0.00 1344(1  O'23) 

Same  as 

15.0 

-258 

0.01463(10‘23) 

RDX  (S.C.)*** 

78.0 

-195 

0.1 629(1  O'23) 

I 

100.0 

-173 

0.2000(1  O'23) 

i 

150.0 

-123 

0.261 1(10'23) 

200.0 

-73 

0.3222(1  O'23) 

218.0 

-55 

0.3442(1  O'23) 

250.0 

-23 

0.3833(  1 0"23) 

See  b 

ote  1 

293.0 

20 

0.4359(1  O'23) 

300.0 

27 

0.4444(1  O'23) 

348.0 

75 

0.503 1(1  O'23) 

350.0 

77 

0.5056(10'23) 

400.0 

127 

0.5667(10'23) 

450.0 

177 

0.6277(1  O'23) 

500.0 

227 

0.6889(10'23) 

525.0 

252 

0.7 194(1  O'23) 

550.0 

277 

0.7500(10'23) 

1 

f 

600.0 

327 

0.7600(1  O'23) 

t 

650.0 

377 

0.7600(1  O'23) 

Same  as 

669.0 

396 

0.7600(1  O'23) 

TATB 

700.0 

427 

0.7600(1  O'23) 

1 

***S.C.  =  Single  Crystal 


Note  1:  CP  [Cal/(Atom  -  °K)  =  [0.2000  +  0.001222  (T  -  100)]  10'23 

=  [0.07780  +  0.001222  T  (°K)]  10'23 
For:  100  <  T  (°K)  <  550 
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Table  2.  FOX-7  Cp  Prediction  Via  the  Proposed  N/G  Cp 


T 

T 

CP 

(N/G) 

NAPMW 

CP 

FOX-7 

CP 

FOX-7 

CP 

FOX-7 

°K 

°C 

Cal. 

Atoms 

Cal. 

Joules 

Joules 

Atom  -  °K 

MW 

MW  -  °K 

MW  -°K 

Gram  -  °K 

0 

-273 

0.0 

84.315  (10+23) 

0.0 

0.0 

0.000 

100 

-173 

0.2000  (10‘23) 

84.315  (10+23) 

16.863 

70.555 

0.476 

150 

-123 

0.2611  (10‘23) 

84.315  (10+23) 

22.015 

92.109 

0.622 

200 

-73 

0.3222  (10‘23) 

84.315  (10+23) 

27.166 

113.664 

0.768 

250 

-23 

0.3833  (10‘23) 

84.315  (10+23) 

32.318 

135.218 

0.913 

300 

27 

0.4444  (10‘23) 

84.315  (10+23) 

37.470 

156.773 

1.059 

350 

77 

0.5056  (10‘23) 

84.315  (10+23) 

42.630 

178.363 

1.204 

400 

127 

0.5667  (10‘23) 

84.315  (10+23) 

47.781 

199.917 

1.350 

450 

177 

0.6277  (10‘23) 

84.315  (10+23) 

52.925 

221.436 

1.495 

500 

227 

0.6889  (10‘23) 

84.315  (10+23) 

58.085 

243.026 

1.641 

550 

277 

0.7500  (10‘23) 

84.315  (10+23) 

63.236 

264.580 

1.787 

600 

323 

0.7600  (10‘23) 

84.315  (10+23) 

64.079 

268.108 

1.810 

650 

377 

0.7600  (10‘23) 

84.315  (10+23) 

64.079 

268.108 

1.810 

MW  (FOX-7)  =  148.086  Grams/Mole 
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Table  3.  FOX-7  Computed  Specific  Heat  Results  for  FOX-7  From  Reference  8,  (Fig.  6) 


T 

T 

CP* 

(Fig. 6) 

CP 

CP 

CP** 

°K 

°C 

Joules 

Joules 

Cal. 

Cal. 

Cell  Mole -°K 

Mole  -  °K 

Mole  -  °K 

Atom  -  °K 

0 

-273 

0.00 

0.00 

0.0000 

0.0000(1  O'23) 

50 

-223 

180.00 

45.00 

10.7553 

0.1276  (10‘23) 

100 

-173 

300.00 

75.00 

17.9254 

0.2126  (10‘23) 

150 

-123 

380.00 

95.00 

22.7055 

0.2693  (10‘23) 

200 

-73 

475.00 

118.75 

28.3819 

0.3366  (10‘23) 

250 

-23 

550.00 

137.50 

32.8633 

0.3898  (10‘23) 

300 

(RT) 

27 

640.00 

160.00 

38.2409 

0.4535  (10‘23) 

350 

77 

705.00 

176.25 

42.1248 

0.4996(1  O'23) 

400 

(EOD) 

127 

770.00 

192.50 

46.0086 

0.5457  (10‘23) 

*  The  Cp  from  Fig.  6  of  Reference  [8]  is  for  4  moles  of  FOX-7. 
FOX-7  — >  C2H4N4O4  — ►  q  =  14.0  Atoms/Molecules 
NAv  =  Avogadro’s  number  =  6.02252  (10+23)  =  Molecules/Mole 
NAPMW  =  q  Nav  =  84.3 15  (10+23)  =  Atoms/Mole 
**  CP  [Cal/  (Atom  -  °K)]  =  CP  [Cal/  (Mole  -  °K)]  /  (NAPMW) 

RT  =  Room  Temperature 
EOD  =  End  of  Data 
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Table  4.  FOX-7  Cp  (T)  Prediction  Comparison 


T 

T 

Cp, 

CP2 

ACP 

A>o) 

°K 

°C 

Joules 

Joules 

Joules 

Percentage 

Difference 

Mole  -  °K 

Mole  -  °K 

Mole  -  °K 

Ref.  [8], 
(Table  3) 

Present 
Report 
(Table  2) 

Cp2  -  Cp, 

% 

0 

.273 

0.000 

0.000 

0.000 

0.00 

100 

-173 

75.000 

70.555 

-4.445 

-5.93 

150 

-123 

95.000 

92.109 

-2.891 

-3.04 

200 

-73 

118.750 

113.664 

-5.086 

-4.28 

250 

-23 

137.500 

135.218 

-2.282 

-1.66 

300  (RT) 

27 

160.00 

156.773 

-3.227 

-2.02 

350 

77 

176.250 

178.363 

+2.113 

+  1.20 

400  (EOD) 

127 

192.500 

199.917 

+7.417 

+3.85 

ACp 

~gT 


(100) 


(  Cpj  -  Cpi 
V  C.P,  , 


x  100.00  =  Percent  Difference 
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Absolute  Temperature,  T,  °K 


Figure  1.  FOX- 7  Cp  Prediction  Comparison 


Figure  2.  Comparison  of  the  Proposed  N/C  Cp  and  the  Calculated  Cp  for  FOX-7 

From  Reference  8 
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APPENDIX  A 

COMPUTATION  OF  mAv,  diav  and  MW  for  FOX-7 


APPENDIX  A 


Computation  of  mAv,  diav  and  MW  for  FOX-7 

The  solid  materials  considered  in  this  study  were  chemical  mixtures.  For  these 
mixtures,  the  weighted  average  mass,  niAv,  of  a  single  atom  in  the  material  was  desired. 

First,  it  was  necessary  to  compute  the  mass  of  a  single  atom  for  each  of  the 
elements  contained  in  the  solid.  The  solid  was  composed  of  one  or  more  of  the  following 
elements: 


Carbon,  C;  Hydrogen,  H;  Nitrogen,  N;  Oxygen,  O: 


The  mass  of  a  single  atom  of  these  elements  is: 


AW  gram/(gram-  mole)  grams 
Nav  atoms/(gram-mole)  atom 


[A-l] 


Where: 


AW= 


grams 
gram  -  mole 


O  C 

Nav  =  Avagadros  Number  =  6.02252  x  10 1 23  - 

gram  -  mole 

Table  A-l  lists  AW  and  m  for  each  of  the  elements  in  the  above  list.  Values  of 
Nav  and  AW  are  from  various  chemistry  text  books  and  handbooks. 

To  compute  the  average  weight  (niAv)  of  an  atom  in  the  material,  the  chemical 
formula  or  proportional  chemical  composition  is  required.  Of  course,  the  weight  (m)  of 
each  elemental  atom  must  be  known,  since  mAv  is  just  a  weighted  average  of  the 
elemental  atoms  in  the  material.  The  procedure  is  valid  for  mixtures  of  compounds  as 
well  as  compounds.  See  References  [1  &  3]  for  examples. 

When  mAv  is  computed,  then  the  average  space  between  the  atoms  (diAv)  is  given 
by  the  following  relation: 

diav  =  ( - Y  =  cm  [A-2] 

P 

Computations  of  mAv  and  diAv  for  FOX-7  are  included  in  this  appendix.  The 
Molecular  weight,  MW,  is  also  calculated. 


Computation  of  mAv  for  FOX-7 


Table  A-l.  Mass  of  a  Single  Atom  for  Selected  Elements 


Element 

AW 
Grams 
Gram  -  Mole 

Nav 

Atoms 
Gram  -  Mole 

m 

Grams 

Atom 

Carbon  (C) 

12.011 

6.02252(1023) 

1.9943(1  O'23) 

Hydrogen  (H) 

1.008 

6.02252(1023) 

0.1 674(1  O'23) 

Nitrogen  (N) 

14.008 

6.02252(1023) 

2.3259(10'23) 

Oxygen (0) 

16.00 

6.02252(1023) 

2.6567(10'23) 

FOX-7  chemical  composition:  C2  H4N4  O4  [8,  11] 


Table  A-2.  FOX-7  mAv  Computation 


C2  2  x  1.9943  (1  O'23) 
H4  4x  0.1674  (1  O'23) 
N4  4  x  2.3259  (10'23) 
04  4  x  2.6567  (10'23) 


=  3.9886(1  O'23)  Grams 
=  0.6696(1  O'23)  Grams 
=  9.3036  (10'23)  Grams 
=  10.6268  (1  O'23)  Grams 


14  Atoms  weigh  24.5886  (10'?3)  Grams 


mAV 


24.5886  (10~23)  Grams 
14  Atoms 


1.756329  (10'23) 


Grams 

Atoms 


23 

«  1.7563  ( 10'“  )  Grams  per  average  Atom 


Computation  of  diAv  for  FOX-7 


Reference  [11]  gives  the  following  density  values  for  FOX-7. 
/?o  =  1.65  Grams/CC  (for  small  scale  gap  test) 

po  =  1.878  Grams/CC  (for  crystal  density) 


For  po  =  1.65  Grams/CC 


dl4V3  = 

raw  _  17.563  (10  24) 

Grams 

U1AV 

P 

1.65 

Grams/CC 

= 

10.644242  (10‘24) 

CM3 

dlAV  = 

2.199741  (10‘8) 

CM 

2.1997  A 

►rj 

O 

>-s 

> 

II 

1.878  Grams/CC 

d,.,,3  = 

niAv  _  17.563  (10  24) 

Grams 

UlAV 

P 

1.878 

Grams/CC 

= 

9.351970  (1  O'24) 

CM3 

dlAV  = 

2.106852  (10‘8) 

CM 

2.1069  A 

CM3 


CM3 


These  values  for  d^v  are  much  larger  than  C(l)  -  C(2)  bond  lengths  (1.456  A)  and 

N(12)  -  0(22)  bond  lengths  (1.242  A  )  given  in  Reference  [11].  However,  Figure  6  of 
Reference  [11]  shows  a  layered  structure  of  the  molecules  where  the  minimum  separation 

distance  between  nitrogen  and  oxygen  atoms  is  3. 1 1  A  .  So  the  average  distance  between 

atoms  is  approximately  2. 10  A  as  shown  in  the  above  computations. 


Computation  of  MW  for  FOX-7 


Chemical  formula  =  C2  H4  N4  O4  =  Ci  Hj  Nk  Oi 

Table  A-3.  MW  Calculation  for  C2  H4  N4  04 


Element 

AW 

n 

n  x  AW 

grams 

grams 

gram  -  mole 

gram  -  mole 

Carbon  (C) 

12.011 

i  =  2 

24.022 

Grams 

Hydrogen  (H) 

1.008 

J  =  4 

4.032 

Grams 

Nitrogen  (N) 

14.008 

k=  4 

56.032 

Grams 

Oxygen  (O) 

16.000 

1  =  4 

64.000 

Grams 

MW  =  Zn  x  AW 

= 

148.086 

grams 

mole 

So  MW  =  148.086  Grams  «  148  Grams 


NAPMW 


MW  =  148.086  Grams 

mAv  1 .756329(1 0  23 )  ( Grams" 

^  Atom  , 

=  84.315638  10  Atomspermole 


This  checks  well  with: 

NAPMW  =  q  Nav  =  (14.0)  (6.02252  *  10+23) 

=  84.315280  *  10+23 


Atoms 


Where: 


q  = 


Zn  =  i  +  j  +  k  +  l  =  2  +  4  +  4  +  4  =  14 


Atoms 

Molecule 


NAv  =  Avogadro’s  number  -  No.  of  molecules  per  MW. 
(Molecules/MW) 
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Appendix  B 

Specific  Heat  (Cp)  Conversion  Relations 


Experimental  Specific  heat  data  are  usually  given  in  calories  (cal)  or  Joules  (J)  per 
molecular  weight  (MW)  per  °  centigrade  (°C)  or  °  Kelvin  (°K).  That  is: 

Cp  =  Cp  (cal.  per  MW  per  °  C  or  °K) 
or 

Cp  =  Cp  (J  per  MW  per  °C  or  °K). 

Occasionally,  Cp  is  presented  as  cal.  or  J  per  gram  °C  or  °K. 

Cp  =  Cp  (cal  per  gram  per  °C  or  °K) 


or 


Cp  =  Cp  (J  per  gram  per  °C  or  °K) 

Note  that: 

Cp  (cal.  or  J  per  gram  per  °C  of  °K)  =  Cp  (cal,  or  J  per  MW  per  °C  or  °K)/MW. 

Since  interest  is  in  Cp  (cal  or  J  per  atom  per  °C  or  °K)  then,  Cp  (cal.  or  J  per  atom  per 
°C  or  °K  =  mAv  Cp  (cal.  or  J  per  gram  per  °C  or  °K). 

See  Appendix  A  for  procedures  to  determine  mAv.  Note  that  1°  C  =  1°  K,  and  that  the 
conversion  factor  between  Joules  and  calories  is  4.184  so  that: 

Cp  (J  per  atom  per  °K)  =  4. 184  *  Cp  (cal  per  atom  per  °K). 

Also,  since  1  Joule  =  107  ergs  =  107  grams  (cm/sec)  2 

then, 


Cp  grams  (cm2/sec2)  per  atom  per  °K  =  107.  (J  per  atom  per  °K). 

The  above  relations  were  employed  in  Reference  2,  and  the  present  report.  However, 
Cp  can  also  be  given  in  terms  of  a  velocity  squared  per  °K  as  follows: 

2  2  7 

Cp(cm  /sec"  per  °K)  =10  x  Cp  (J  per  gram  per  °K).  In  these  units  for  Cp,  A(v.e.)xR  = 
f  CpdT  =  (cm/sec)2,  so  Upcri  =  J A(v.c.)tr  =  and  Upcr:  =  J 2A(v.c)tr  =  V2  Upcri  . 

JTexp  v  v 

This  eliminates  having  to  compute  an  average  mass  (mAv)  and  the  Cp  in  calories 
or  Joules  per  average  atom. 
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